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Experiments have been performed using a spherical superconducting microwave resonator that 
simulates the geometric structure of the Cso fullerene molecule. The objective was to study with very 
high resolution the exceptional spectral properties emerging from the symmetries of the icosahedral 
structure of the carbon lattice. In particular, the number of zero modes has been determined to test 
the predictions of the Atiyah-Singer index theorem, which relates it to the topology of the curved 
carbon lattice. This is, to the best of our knowledge, the first experimental verification of the index 
theorem. 

PACS numbers: 05.45.Mt,41.20.Jb,71.20.-b,71.20.Tx,73.22.-f 


Introduction .— The spectrum of graphene, a mono- 
layer of carbon (C) atoms arranged on a hexagonal lat¬ 
tice, has been the focus of extensive theoretical [l|,y and 
experimental studies Q. Its universal properties were 
often also investigated experimentally in analog systems, 
so-called ’artifici al g raphene’ Q, e.g., in our group in pho¬ 


tonic crystals [5l-ll0j. Moreover, theoretically much atten¬ 
tion has been devoted to curved graphene structures like 
fullerene molecules [lll - [l6l| and the connection between 
their spatial symmetries and electronic properties. Here, 
the most famous example is the Ceo molecule. It con¬ 
sists of 60 carbon atoms at the vertices of a truncated 
icosahedron and has the shape of a soccer ball. Con¬ 
cerning the spectral properties of fullerenes the number 
of near-zero modes, i.e., of electronic states with excita¬ 
tion energies close to zero, have been of particular in¬ 
terest since they determine the electrical conductivity. 
In [13, [3 an index theorem has been derived that allows 
the computation of the number of such near-zero modes 
from the topology of the surface. It was deduced from 
the renowned Atiyah-Singer index theorem l^- which 
states that the analytic index of an elliptic differential op¬ 
erator on a compact manifold equals the topological one, 
in other words, that there is a connection between the 
number of zero modes of the operator and the topology 
of the manifold on which it is defined. The aim of the 
high-resolution experiments presented in this letter was 
to test these predictions in experiments with a supercon¬ 
ducting microwave resonator of the same topology as the 
Ceo molecule. 


First we briefly review the salient features of graphene 
and fullerenes and outline the derivation of the index 
theorem from 33 for deformed graphene sheets. We 
then describe the experimental setup and compare the 
results of the measurements to the predictions from the 
index theorem and to tight-binding model (TBM) calcu¬ 
lations. These allow us to study the approach to the ther¬ 
modynamic limit of an infinite number of carbon atoms. 


Graphene, fullerenes and the Atiyah-Singer theorem .— 
The honeycomb structure of graphene is formed by 
two interpenetrating triangular sublattices. As a con¬ 
sequence, at half filling the Fermi surface in graphene 
reduces to two independent points in the first Brillouin 
zone, the so-called ’Dirac points’, denoted by K_(. and 
K_ [l|,l4| that are conical intersections of the valence and 
the conduction band. Low energy excitations within the 
cone regions around K± have a linear dispersion with 
a slope given by the Fermi velocity vp- On an infinte 
graphene sheet they are therefore described by a Dirac 
Hamiltonian for massless spin-1/2 quasiparticles consist¬ 
ing of partner Hamiltonians 

H± = ±VF <J°‘qa (1) 


which describe excitations with momentum q = iQxiQy) 
in each of the two Dirac cones around K±. The Pauli 
matrices (t“ with a = x,y act on the two sublattice com¬ 
ponents of the excitations, combined in two-dimensional 
spinors and hence referred to as quasi-spin. Both cones 
together then yield a four-component Dirac equation [ 3 ] . 

Fullerene molecules can be constructed by introduc¬ 
ing positive curvature into an initially flat graphene 
sheet 3 ■ The bending is realized by replacing hexagons 
by pentagons, ensuring at the same time that the lattice 
is not stretched and each C atom keeps three neighbors. 
To determine the number of pentagons ns necessary to 
generate a spherical fullerene molecule with uq hexagons 
one uses the Euler formula 2^ which relates the number 
of vertices V, of edges E, faces F and open ends A^open 
of an arbitrary two-dimensional lattice to the genus g of 
the surface formed by it, via the Euler characteristic 


^ = V-E + F = 2{l-g)- A^open ■ 

For a lattice of pentagons and hexagons, V = (Sns -|- 
6ne)/3, E = (Sns -|- 6n6)/2 and F = uq, this gives 
X = ^ 75 / 6 . Without open ends y must be an even integer 




2 


on a closed orientable surface due to the Gaufi-Bonnet 
theorem 2^. Hence, for a flat graphene sheet with peri¬ 
odic boundary conditions one has 5 = 1 for the torus and 
ns = 0 , while a sphere with 5 = 0 needs = 12 pen¬ 
tagons to avoid open ends. Consequently, fullerenes are 
grown from the Ceo molecule by increasing the number 
of hexagons, i.e., always have twelve pentagons at the 
same relative positions. This also applies to the ther¬ 
modynamic limit, and one expects that their low-energy 
electronic excitations are described by a Dirac equation 
on a sphere. 

To introduce a pentagon into the honycomb lattice, a 
7 r /3 sector is cut out and then the edges are glued to¬ 
gether 13, 3, 12, as illustrated in Fig. [T] Thereby 



FIG. 1: (Color online) Left panel: A two-dimensional 
graphene sheet. The red and the blue dots mark the two 
independent triangular sublattices. In order to form a curved 
sheet which contains one pentagon, a 7 r /3 segment is cut 
out from the sheet. Right panel: The conically deformed 
graphene sheet with the pentagon at the apex. 


a pentagon is created at the apex of the emerging cone. 
Along the seam, two C atoms from the same triangular 
sublattice, e.g., the red ones in Fig. [U are connected. 
This results in a coupling of the Dirac operators as¬ 
sociated with the K± points. Indeed, when the four¬ 
dimensional spinor associated with the Dirac equation of 
the flat graphene sheet is transported around the apex 
by an angle 27r, it is forced to jump at the seam from 
a red site to another red one instead of to a blue one. 
It thus acquires a non-trivial phase, which can be ac¬ 
counted for by introducing a non-Abelian gauge field 
in the Hamiltonian which yields a flux of ( 7 r/ 2 ) when 
integrated along a closed loop around the apex. The 
Pauli matrix thereby couples the K+ and K_ spinor 
components [ij. This description entails the existence 
of a ficitious magnetic monopole inside the surface. In 
the case of fullerenes it is located at the center of the 
spherical molecule, yielding a flux of 1 /8 through each of 
the twelve pentagons. Thus the total magnetic monopole 
charge inside the sphere equals 3/2. In addition to that, 
analogous to the daily rotation of Foucault’s pendulum, 
a deficit angle of 7 r /3 arises when moving a frame along 
a loop around the apex. It is here described by a quasi¬ 


spin connection with circulation — ( 7 r/ 6 )cr^ around 
the apex. 

The coupling of the K± spinor components in the re¬ 
sulting four-dimensional Dirac equation can be removed 
by a rotation, which leads to two independent two- 
dimensional Dirac equations denoted by Z = 1, 2 |18l| . 


= VF o-“e^ ( 5 ^ - iQ^ - iA‘^) ip'' = Eip' 


( 2 ) 


where is the Zweibein in the tangent plane of the 
surface, and A'^ are the components of in the rotated 
basis with circulation ±7r/2 for I = 1 and 2, respectively, 
where A' is now an Abelian gauge field. 

The four-dimensional Dirac equation obtained from 
Eq. ([2]) provides a good description of the low-energy 
excitations of the C molecules l3-15|- It yields the long- 
wavelength excitations of the deformed graphene sheet 
in the vicinity of the Dirac points, and thus also the zero 
modes that we are interested in. For the fullerenes the 
Dirac operators ij) are elliptic and defined on a compact 


surface. Hence the Atiyah-Singer index theorem [I9l-l2 


applies. Ten years after its first formulation a new proof 
was provided based on the heat equation [ 3 ] which was 
later employed for the derivation of an index theorem for 
graphene sheets deformed by pentagons and heptagons 
[171 [l 8 j as briefly reviewed in the following. 

Each Dirac operator in Eq. (ED can be written in terms 
of off-diagonal partner operators P and [3 1 so (ip'r 
contains only diagonal operators PP^ and pIp that have 
the same number of zero modes as Pl and P, respectively. 
Furthermore, the non-zero eigenvalues of PPl and pIp 
are identical. The analytic index of Ip is given by the 
difference of the numbers of zero modes of P and pl 
denoted by v±, respectively, i.e., index( 2 ^ ) = 1 ^+ — v- 
17 . 3 . More importantly, however, this index is related 
to the total flux of the effective gauge field via the Atiyah- 
Singer index theorem 


index( 0 *) = ^ // . 


(3) 


The integral is taken over the compact surface D and 
P' = dAA' are the field strengths associated with the 
now Abelian gauge potentials A'. Stokes’ theorem then 
implies from the closed loops around each apex that 



The Euler formula thus leads to the Atiyah-Singer index 
theorem for fullerenes in the form, index(0*) = ± 3 ( 1 — 5 ). 
In two dimensions either or vanish. Hence, the 
index theorem provides the number of zero modes [3 . 

m. 

The total number of zero modes is the sum of those of 
the subsystems corresponding to I = 1,2. Consequently, 
according to the index theorem, the zero modes of the 
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four-dimensional Dirac operator for spherical fullerenes 
correspond to two triplets. The same result has been 
obtained in a continuum model for the low-energy elec¬ 
tronic states of icosahedral fullerenes 13 We em¬ 
phasize, however, that the eigenvalues of the near-zero 
modes tend to zero, i.e., coincide with the energy at the 
Dirac points, only in the thermodynamic limit of an in¬ 
finity number of C atoms. In a sufficiently large but finite 
fullerene molecule they are expected to lie much closer to 
the Dirac energy than all the other ones. 

Experimental setup and resonance spectra .— Hitherto, 
experiments have been performed with flat, supercon- 
ductin g m icrowave resonators, so-called ’microwave bil¬ 
liards’ [27| to address problems from the fields of quan¬ 
tum chaos 28, l^l and compound nucleus reactions [3nj |. 
In this context, the equivalence of the Helmholtz equa¬ 
tion and the non-relativistic Schrodinger equation of the 
corresponding quantum billiard is exploited which holds 
below a maximum microwave frequency /max = c/(2d) 
with c the velocity of light and d the height of the bil¬ 
liard. Consequently, the eigenvalues of a quantum bil¬ 
liard can be obtained experimentally from the eigenfre- 
qencies of the microwave billiard of corresponding shape. 
Recently, we realized experiments with superconducting 
microwave Dirac billiards and studied universal spectral 
properties of graphene sheets [sij with unprecedented ac¬ 


3 , 132 , 133 . 


curacy 

The aim of the experiments presented here was the 
investigation of the universal spectral properties of the 
fullerene Ceo molecule attributed to its lattice structure 
and to determine the number of zero modes which, ac¬ 
cording to the Atiyah Singer index theorem solely de¬ 
pends on the number of pentagons. For this we use a sys¬ 
tem exhibiting the same topological properties, namely 
a quantum fullerene billiard on a sphere, consisting of a 
network of 60 circular billiards at the positions of the C 
atoms connected by three of the altogether 90 straight 
leads with three adjacent ones. We studied them exper¬ 
imentally by using instead of a planar, superconducting 
microwave (Dirac) billiard a cavity, which is imprinted 
on a sphere. The microwave fullerene billiard displayed 
in Fig. [2] was constructed by milling a total of 60 circu¬ 
lar cavities (vertices) and 90 rectangular channels (edges) 
out of a brass sphere and then closing them with small 
triangular brass plates of 5 mm thickness and 3 mm thick 
rectangular ones, respectively. Before the parts were 
screwed together, they were covered with lead, which is 
superconducting below Tc = 7.2 K. The diameter of the 
sphere of 160 mm was limited by the size of the liquid 
Helium cryostat in which the resonator was cooled down 
to 4.2 K in order to attain superconductivity. The ra¬ 
dius of the circular cavities was 12 mm, the widths of 
the waveguides 14 mm, before lead coating them. Thus 
the cutoff frequency for the hrst propagating mode in the 
latter is // > 10.714 GHz. In total, 8 antennas were at¬ 
tached to the triangular plates. Two, covered with red 



FIG. 2: (Color online) Lead plated fullerene billiard used in 
the experiments. In the left part the small plates that cover 
the circular cavities and the rectangular channels were re¬ 
moved. The red caps protect the antenna ports. The billiard 
is superconducting below Tc=7.2 K. 


caps, are visible in Fig. [T] The height of the resonator 
was 3 mm corresponding to /max = 50 GHz. 

For the measurement of the transmission spectrum 
shown in the upper panel of Fig. [3] microwave power was 
coupled into the microwave billiard via antenna a and 
the output signal was received at antenna b, with a and 
b denoting two of the 8 antennas. A vectorial network 
analyzer determined the relative phase and amplitude of 
the output and input signals, thus yielding the scattering 
matrix element Sha- The smallest resonance frequency 
equals / = 8.254 GHz, so we show the spectrum from 
8-40 GHz. Due to the high-quality factor Q > 10^ of 
the resonator, all resonances could be resolved in that 
frequency range. We concentrate our discussion here on 
the region between 8.254 and 18.801 GHz. The spec¬ 
trum exhibits three distinct bands, containing 60, 210 
and 90 resonances, respectively, in the frequency intervals 
[8.254, 8.779] GHz, [11.492, 16.657] GHz and [18.312, 
18.801] GHz. They are located around the eigenfrequen- 
cies /i ~ 8.4 GHz, /a ~ 13.5 GHz and /a ~ 18.5 GHz 
of the first three quasibound states in an open circular 
billiard of the same size as the cavities in the resonator 
with openings at the positions of the waveguides. The 
modes excited inside the circular cavities resemble within 
a given band the corresponding mode in the open circular 
billiard, and are described by Jq, Ji and Ja Bessel func¬ 
tions in the first, second and third band, respectively. 
In the latter two cases they are twofold degenerate due 
to the mirror symmetry. Note, that the circular cavities 
exhibit no threefold symmetry, because each of them is 
part of one pentagon and two hexagons and the internal 
angles differ. 

The first band is located well below the cutoff fre¬ 
quency of the waveguide. Consequently, the electric field 
modes excited inside the circular cavities are only weakly 
coupled to those in the neighboring ones. The resonance 
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FIG. 3: Transmission spectrum of the fullerene billiard (upper 
panel) up to 40 GHz. The first band ranges from 8.254 to 
8.779 GHz. It contains 60 resonances, that are separated 
into 15 groups with the number of resonances indicated in 
the middle panel, which shows a zoom into it. The Dirac 
frequency fn = 8.504 GHz is marked by an arrow. The two 
triplets of interest (the zero modes) are clearly resolved and 
shown in the lower panel. 


frequencies in the second band are above the cutoff fre¬ 
quency. Accordingly, the modes in the cavities are cou¬ 
pled via the modes inside the waveguides, and thus mim- 
ick a situation where the C atoms are coupled to the 
neighboring ones via an extra atom, thus explaining the 
number of resonances in this band. The third band is 
still below the frequency > 20.143 GHz of the sec¬ 
ond propagating mode in the waveguides. As a result, 
the number of possible mode configurations is restricted 
due to the symmetry properties of the modes excited in¬ 
side the cavities, that prefereably couple to the second 
excited mode inside the waveguide. Above 20.232 GHz, 
i.e., beyond /^, several bands are intertwined. 

In summary, only the first band can be used to model 
the situation in the fullerene Geo molecule. The middle 
panel of Fig.[3]shows a magnification of it. Fifteen groups 
of nearly degenerate resonances are clearly visible. The 
number of resonances identified in each of them is indi¬ 
cated and coincides with the degrees of degeneracy pre¬ 
dicted on the basis of group theoretical considerations for 
the eigenfrequencies because of its truncated icosahedral 
structure [l^, IsJ • In the group with degeneracy degree 
9, in fact, the energy values of 5 and 4 degenerate eigen¬ 
frequencies, respectively, are accidentally the same. We 
emphasize, that we were only able to identify all 60 res¬ 
onance frequencies because the degeneracies were lifted. 
The reason is that the symmetries of the resonator struc¬ 
ture were slightly perturbed due to the presence of the an¬ 
tennas and unavoidable marginal inhomogeneities in the 


lead coating. The influence of the former turned out to 
be negligible for sufficiently short antennas. The effect of 
the latter on the size of the splittings of the nearly degen¬ 
erate resonance frequencies was tested by smoothing the 
surface of the resonator which indeed induced a reduction 
of the splittings in each group. The pair of triplets visible 
in the middle panel of Fig. [3]and shown in a further mag¬ 
nification of the spectrum in the lowest panel, is closest 
to the Dirac frequency at fo = 8.504 GHz which was 
determined as described below. The zoom demonstrates 
the high resolution necessary to resolve the 6 resonances. 
These are the 6 modes conjectured by the Atiyah-Singer 
index theorem that we were looking for. As is discussed 
next they are corroborated also by TBM calculations. 

Tight-binding model description of the spectra .— The 
eigenvalues of the Geo molecule have been computed pre¬ 
viously using the TBM 12|, however, a stringent test 
of its applicability was missing. Given the experimental 
results on the eigenfrequencies in the first band of the 
fullerene resonator we are now in a position to check 
the validity of the TBM in detail. As stated above, 
the modes excited in the 60 cavities are weakly coupled, 
which is an essential prerequisit for the applicability of 
the TBM. Detailed calculations showed a quantitative 
agreement between the computed and the measured fre¬ 
quencies only when including next-nearest, and second 
and third-nearest neighbor couplings with strengths G, t 2 
and ta, respectively. This yielded for the frequency of 
the isolated cavities /o = 8.515 GHz and the coupling 
parameters ti = —0.0929 GHz, t 2 = 0.0035 GHz and 
ta = 0.0005 GHz. The eigenvalues deduced from the 
TBM appear as 15 groups of degenerate eigenvalues with 
the same multiplicities as the resonances in the spectrum 
shown in the middle panel of Fig. [3l An even better 
agreement was achieved by taking into account the fact 
that, due to the inhomogeneities in the lead coating, the 
radii of the cavities are slightly different. In order to esti¬ 
mate the deviations thus induced in /o, we used the fact, 
that /o is given by the first zero of the Jq Bessel func¬ 
tion, Jo{kR) = 0, for a circular cavity of radius R. We 
inserted for each of the 60 cavities the measured radius 
and replaced /o by the individual values. Thereby, the 
degeneracies were removed. In panel a) of Fig. |4] we com¬ 
pare the resonance density p{f) = ^ ■ 6{f — fi) obtained 
for the frequencies fi in the first band (black full line) 
with the TBM result (red dashed line). For display pur¬ 
poses we have replaced the S functions by Lorentzians of 
finite width of F = 2 MHz. The good agreement reassures 
the applicability of the TBM and thus justifies its use for 
further numerical studies with larger fullerene molecules 
with the parameters determined from the experiment. 

Panel b) of Fig.|4]shows a comparison between the reso¬ 
nance densities of the Ceoj G240, C540 and C720 molecules 
in ascending order. Here, we used Lorentzians of width 
F = 5 MHz for all cases. As stated above, all molecules 
contain the same number of pentagons whereas the num- 
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FIG. 4: (Color online) Panel a): Resonance density deter¬ 
mined from the resonance frequencies of the first band (mid¬ 
dle panel in Fig. 0 of the fnllerene billiard (black solid line) 
compared to the TBM calculations (red dashed line). Panel 
b): Resonance densities compnted within the TBM for Ceo, 
C 240 , C 540 and C 720 (lowermost to upmost cnrve). Panel c): 
Experimentall y d etermined resonance density of a rectangnlar 
Dirac billiard (lOl. 1^. 


ber of hexagons increases. The resonance densities should 
thus resemble more and more that of a graphene sheet. 
They exhibit a minimum bounded by two increasingly 
sharp peaks, that evolve into van Hove singularities in 
the limit of an infinite number of atoms sH], 36| . A com¬ 
parison of the resonance densities with that of a rectan¬ 
gular graphene sheet with periodic boundary conditions 
(panel c) of Fig. 0), which has been obtained in mea¬ 
surements with a microwave Dirac billiard [l^ . shows 
that they resemble for large fullerenes. In contradistinc¬ 
tion to the latter, however, the resonance densities of the 
fullerenes all exhibit a peak of similar size located at the 
minimum which is due to the two triplets of zero modes. 
This remains true in the limit of an infinite number of 
atoms [IMl and is thus a distinct feature of the spa¬ 
tial curvature and topology. According to the Atiyah- 
Singer index theorem a plane graphene sheet with pe¬ 
riodic boundary conditions should not exhibit any zero 
modes. This is observed in panel c). Zero modes are ex¬ 
pected only, if part of the graphene sheet is terminated 
with zigzag edges [3l|. The associated states are called 
edge states, because their wave functions vanish every¬ 
where except at these edges 1^ . We have also 

computed the wave functions of the fullerene molecules 
under consideration using the TBM and found, that those 
of the 6 zero modes are localized at the pentagons, which 
may be considered to be equivalent to zigzag edges within 
the hexagon network. 

The central frequencies of the two triplets, marked by 
squares for the one close to the Dirac frequency (dotted 


line) and by a circle for the one further away are displayed 
in Fig. [5] as a function of the number n of C atoms. As is 
clearly visible, the distance between the triplets decreases 
with increasing size of the fullerene molecule and both 
approach the Dirac frequency. This behavior is well fitted 
by a function /(n) = fo + a/n}’ yielding the parameter 
values given in the caption of Fig. [Sj and for the Dirac 
frequency finally a value of fo = 8.504 GHz. Note that 
this is essentially the only way to determine the Dirac 
frequency of a Ceo molecule. 



Number n of C Atom.s 


FIG. 5: (Golor online) Calculated frequencies of the zero 
modes of the fullerenes Ceo, Ciso, C 240 , C 540 and G 720 vs. 
the number n of C atoms. All pairs of triplets (squares and 
circles) are located slightly above the Dirac frequency (dot¬ 
ted line). The experimental values are marked by a red filled 
circle and square. The dashed lines correspond to fits of the 
function /(n) = fo + a/n^ to the data points yielding for the 
Dirac frequency fo = 8.504 GHz, a = 1.188 GHz, b = 0.661 
for the lower zero modes and a = 1.211 GHz, b = 0.869 for 
the upper ones. 


Conclusions .— The lowest 60 eigenvalues of a Ceo 
fullerene were determined in high-precision experiments 
using a superconducting microwave billiard of corre¬ 
sponding shape. They appear in 15 groups of nearly de¬ 
generate ones, where the multiplicity coincides with that 
determined based on the group theory of the truncated 
icosahedral structure of Ceo- We have demonstrated in 
TBM calculations for spherical fullerene molecules of in¬ 
creasing size that the two triplets of resonances detected 
close to the Dirac frequency correspond to the triplets of 
zero modes predicted by the Atiyah-Singer index theo¬ 
rem, and thus provided to the best of our knowledge the 
first experimental test of it. The exact value of the Dirac 
frequency was obtained as the asymptotic value attained 
by the frequencies of the triplets in the limit of an infinite 
number of atoms. 

This work was supported by the DFG within the Col¬ 
laborative Research Center 634. 
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